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ABSTRACT 

The flat spectrum radio quasar (FSRQ) PKS 2123—463 was associated in the First 
-Ferrai-LAT source catalog with the 7 -ray source 1FGL J2126. 1—4603, but when con- 
sidering the full first two years of Fermi observations, no 7 -ray source at a position 
consistent with this FSRQ was detected, and thus PKS 2123—463 was not reported in 
the Second Fermi - LAT source catalog. On 2011 December 14 a 7 -ray source position- 
ally consistent with PKS 2123—463 was detected in flaring activity by Ferrm-LAT. 
This activity triggered radio-to- X-ray observations by the Swift, GROND, ATCA, Ce- 
duna, and KAT-7 observatories. Results of the localization of the 7 -ray source over 41 
months of Fermi - LAT operation are reported here in conjunction with the results of 
the analysis of radio, optical, UV and X-ray data collected soon after the 7 -ray flare. 
The strict spatial association with the lower energy counterpart together with a simul- 
taneous increase of the activity in optical, UV, X-ray and 7 -ray bands led to a firm 
identification of the 7 -ray source with PKS 2123—463. A new photometric redshift 
has been estimated as z= 1.46 ± 0.05 using GROND and Swift / UVOT observations, 
in rough agreement with the disputed spectroscopic redshift of z = 1.67. We fit the 
broadband spectral energy distribution with a synchrotron/external Compton model. 
We find that a thermal disk component is necessary to explain the optical/UV emis- 
sion detected by Swift/ UVOT. This disk has a luminosity of ~ 1.8 x 10 46 erg s _1 , and 
a fit to the disk emission assuming a Schwarzschild (i.e., nonrotating) black hole gives 
a mass of ~ 2 x 10 9 M 0 . This is the first black hole mass estimate for this source. 

Key words: galaxies: active - galaxies: nuclei - quasars: general - quasars: individual 
(PKS 2123-463) - gamma rays 


1 INTRODUCTION 

Blazars constitute the most extreme subclass of Active 
Galactic Nuclei (AGN), characterized by the emission of 

* E-mail: filippo.dammando@fisica.unipg.it 


strong non-thermal radiation across the entire electromag- 
netic spectrum and in particular intense 7 -ray emission 
above 100 MeV. The typical observational properties of 
blazars include irregular, rapid, high-amplitude variability, 
radio-core dominance, apparent super-luminal motion, a flat 
radio spectrum, and high and variable polarization at radio 
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and optical frequencies. These features are interpreted as 
resulting from the emission of high-energy particles acceler- 
ated within a relativistic jet closely aligned with our line of 
sight and launched in the vicinity of the supermassive black 
hole harboured by the activ e galaxy felandford fe Rees I 
ll978l;IUrrv fc Padovani1ll995ll . 

Since the advent of the Energetic Gamma Ray Exper- 
iment Telescope (EGRET) on the Compton Gamma-Ray 
Observatory (CGRO), blazars were known to dominate the 
extragalactic high-energy sky. However, EGRET did not 
pinpoint the location of many sources with sufficient pre- 
cision to enable astronomers to associate them with known 
objects, leaving the legacy of a large fraction of unidentified 
sources in 7 rays. The point-spread function and sensitivity 
of the Large Area Telescope (LAT) on board Fermi pro- 
vides an unprecedented angular resolution at high energies 
for localizing a large number of newly found 7 -ray emitting 
sources. Correlated variability observed at different frequen- 
cies can give important information for the identification of 
a 7 -ray source with its low-energy counterpart. 

On 2011 December 14 a 7 -ray flare from a source 
positionally cons istent w ith PKS 2 123— 4 63 w as detected 
by Fermi - LAT (jOrienti fe D’Ammandd l201llL trigg erini 


GR OND and Swift follow-up observations ( Rau et al 
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iD’Ammando et al] [20111 that confirmed contemporaneous 
activity in the optical/UV as well as marginally in X-rays. 

PKS 2123—463 is a bright radio quasar with a lumi- 
nosity at 1.4 GHz Li.4ghz= (1.5±0.2)x10 28 W/Hz (assum- 
ing the redshift estimated in this paper, z = 1.46±0. 
see Sect. On the basis of its spectral index a r ~0.4 
(S(i') ex u~ ar ) between 408 MHz and 4.8 GHz, it was in- 
cluded in the CRATES catalog of flat spectrum objects 
dHealev et al. 11200 7(1. ATCA observations performed almost 
simultaneously at 4.8, 8.6 and 20 GHz during the Australia 
Telescope 20 GHz (AT20G) survey indicated a flattening of 
the radio spectrum at high frequencies to a spectral index of 
a r ~ 0.2. Polarized emission has been detected only at 4.8 
GHz where the pola rizat ion is abou t 4% of the total intensity 
flux density (iMassardi et al.l l2008l l . A Chandra observation 
of PKS 2123—463 in March 2004 has show n the presence 
of a jet-like extended structure in X-rays ([Marshall et al.l 
1201 llf . The redshift of PKS 2123—463 was reported to be 
z = 1.67 (Sav age fe Wrightlll98ll ') based on an objective- 
prism spectrum, subsequently questioned bv I Jackson et al] 
(l2002ll because two possible redshifts (0.48 and 1.67) were 
given from observation of two lines in the spectrum, and 
the motivation for the exclusion of the smaller redshift was 
not provided. This obj ect was a member of the pr e-Fermi 
launch Roma-BZCAT dMassaro et al.|[2009l 'l catalog listing 
candidate 7 -ray blazar s detectab le by Fer mir- LAT but not in 
the CGRaBS catalog dHealev et al.lf2008ii . 

The flat spectrum radio quasar (FSRQ) PKS 
2123- 463 ( R.A.:21h26m30. 7042s, Dec.:-46d05m47.892s, 
J2000; iFev et al.l 1200^4 was associated in the First Fermi- 


(1FGL, 2008 August 4 - 2009 
with the 7 -ray source 1FGL 


LAT s ource catalog ( it t 
July 4; lAbdo et al] l2010al ~) 

J2126. 1—4603, while no association with a 7 -ray source was 


1 The corresponding luminosity distance for z = 1.46T0.05 is d £, 

= 10.6T0.5 Gpc and 1 arcsec corresponds to a projected distance 
of 8.53T0.02 kpc. 


present in the Second Eerrm- LAT source catalog (2FGL, 
2008 August 4 - 2010 August 4: lNolan et al.|[2012l l. although 
a 7 -ray source, 2FGL J2125. 0—4632 with a radius of the 95% 
source location confidence region of 0.17°, at 0.52° from the 
radio position of PKS 2123—463 was reported. Taking into 
consideration the high variability of blazars, the flux of PKS 
2123—463 could have decreased in the second year of Fermi- 
LAT operation. In addition, in the construction of the 2FGL 
catalog, PKS 2123—463 was split int o more t han one c an- 
didate source seed (see Section 4.2 of lNolan et al.ll2012l . for 
details). These two factors may have led to the lack of its 
association with a 7 -ray source. 

In this paper we present the localization over 41 months 
of Eerrm-LAT data of a 7 -ray source associated with the 
FSRQ PKS 2123—463. The correlated variability observed 
in optical, UV, X-ray and 7 rays confirms the identification. 
In addition a new estimation of the redshift of the source by 
means of the fit of simultaneous GROND and Swift/WOT 
data collected soon after the 7 -ray flare is presented. The pa- 
per is organized as follows: in Section 2 we report the LAT 
data analyses and results. In Sections 3 and 4 we report the 
results of the Swift and GROND data analysis, respectively. 
Radio data collected by the ATCA, Ceduna and KAT-7 tele- 
scopes are presented in Section 5. Section 6 presents an esti- 
mation of the redshift of the source. In Section 7 we discuss 
the modeling of the overall SED and draw our conclusions. 
Throughout the paper a A-CDM cosmology with Ho = 71 
km s ^ 1 Mpc~ x , Ha = 0.73, and = 0.27 is adopted. 


2 FERMI - LAT DATA: SELECTION AND 

ANALYSIS 

The Eerrm-LAT is a 7 -ray telescope operating from 20 MeV 
to > 300 GeV. It has a large peak effective area (~ 8000 cm 2 
for 1 GeV photons), an energy resolution of typically ~ 10%, 
and a field of view of about 2.4 sr with an angular resolution 
( 68 % containment angle) better than 1 ° for energies above 
1 GeV. Further deta ils about the Fermi-LAT are given by 
lAtwood et al.l d2009h . 

The Eerrm-LAT data reported in this paper were col- 
lected over the first 41 months of Fermi operation, from 2008 
August 4 (MJD 54682) to 2012 January 4 (MJD 55930). 
During this time the LAT instrument operated almost en- 
tirely in survey mode. The spectral analysis was performed 
with the ScienceTools software package version v9r23pl. 
The Eerrm-LAT data were extracted from a circular Region 
of Interest (Rol) with a 15° radius centred at the radio loca- 
tion of PKS 2123—463. Only events belonging to the ‘Source’ 
class were used. The time intervals when the rocking angle 
of the LAT was greater than 52° were rejected. In addition, 
a cut on the zenith angle (< 100 °) was also applied to reduce 
contamination from the Earth limb 7 rays, which are pro- 
duced by cosmic rays interacting with the upper atmosphere. 
The spectral analyses (from which we derived spectral fits 
and photon fluxes) were performed with the post-launch in- 
strument response functions (IRFs) P7S0URCE_V6 using an 
unbinned maximum likelihood method implemented in the 
Science tool gtlike. 

The background model used to extract the 7 -ray sig- 
nal includes a Galactic diffuse emission component and an 
isotropic component. The model that we adopted for the 


PKS 2123—463: a confirmed 'y-ray blazar at high redshift 3 


2009 2010 90 1 2012 


3C 

CO 

E 

u 

„ 20 

c 

o 

o 

JI 

a. 

° 10 


0 


Figure 1. 2-month integrated flux (E > 100 MeV) light curve of 
PKS 2123—463 obtained from 2008 August 4 to 2012 January 4 
(MJD 54682—55930). The last bin covers 1 month of observation. 
Arrows refer to 2-er upper limits on the source flux. Upper limits 
are computed when TS < 10. 
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Figure 2. Integrated flux (E > 100 MeV) light curve of PKS 
2123—463 obtained from 2011 December 4 to 2012 January 4 
(MJD 55899-55930) with 2-day or 4-day time bins. Arrows refer 
to 2 -it upper limits on the source flux. Upper limits are computed 
when TS < 10. 


Galactic component, the same as used for the 2FGL catalog, 
is given by the file gal_2yearp7v6_v0.fits, and the isotropic 
component, which is the sum of the extragalactic diffuse 
emission and the residual charged particle background, is 
parametrized by the file iso_p7v6source.tx10. The normal- 
izations of both components in the background model were 
allowed to vary freely during the spectral point fitting. 

We examine the significance of the 7 -ray signal from 


2 http: //fermi. gsfc.nasa.gov/ssc/data/access/lat/Background 

Models.html 


the sources using the Test Statistic (TS) based on the 
likelihood ratio test. The Test Statistic, defined as TS = 
2Alog(likelihood) between models with and without the 
source, is a measure of the probability of having a 7 -ray 
source at the localization specified, which compares models 
whose parameters have been adjus ted to maximize th e like- 
lihood of the data given the model dMattox et al.l[l99^ ). The 
source model used in gtlike includes all the point sources 
from the 2FGL that fall within 20° of PKS 2123—463. 
The spectra of those sources were parametrized by power- 
law functions, except for 2FGL J2056.2— 4715 for which we 
used a log-parabola, and 2FGL J2124.6— 3357 and 2FGL 
J2241.7— 5236 for which we used an exponentially cut-off 
power-law in their spectral modeling as in the 2FGL cat- 
alog. We removed from the model the sources having TS 
< 25 and/or fluxes (0.1 100 GeV) below l.OxlO -8 photons 
cnfl 1 s _1 over 41 months and repeated the fit. We tested 
whether two distinct 7 -ray sources (one at the radio posi- 
tion of PKS 2123—463 and one at the 7 -ray position of 2FGL 
J2125. 0—4632) are detected simultaneously by Fermi with 
TS ^ 25 over 41 months of observations. In this scenario, the 
fit yields a TS of 553 for PKS 2123—463 and a TS of 14 for 
2FGL J2125. 0—4632. We thus conclude that PKS 2123—463 
and 2FGL J2125. 0—4632 are the same source and we main- 
tain only PKS 2123—463 in the model. Thus a final fitting 
procedure has been performed with all sources within 10 ° of 
PKS 2123—463 included with the normalization factors and 
the photon indices left as free parameters. For the sources 
located between 10° and 15° we kept the normalization and 
the photon index fixed to the values obtained in the previ- 
ous fitting procedure. The Rol model includes also sources 
falling between 15° and 20° from the target source, which 
can contribute to the total counts observed in the Rol due 
to the energy-dependent size of the point spread function of 
the instrument. For these additional sources, normalizations 
and indices were fixed to the values of the 2FGL catalog dur- 
ing all steps of the fitting procedure. 

The 7 -ray point source localization using the gtf indsrc 
tool over the photons extracted during the period 2008 Au- 
gust 4 - 2012 January 4 results in R.A. = 321.609°, Dec. = 
-46.076° (J2000), with a 95% error circle radius of 0.047°, 
at an angular separation of 0.024° from the radio position of 
PKS 2123-463 (R.A. = 321.628°, Dec. = -46.097°, J2000). 
This implies a strict spatial association with PKS 2123—463. 
The fit with a power-law (PL) model, dN/dE oc (E/Eo)^ r , 
to the data integrated over 41 months of Fermi operation 
in the 0.1-100 GeV energy range results in a TS = 557, 
with an integrated average flux of (4.76 ± 0.35) xl0 ~ 8 
photons cm ' 2 * s _1 , and a photon index T = 2.45 ± 0.05 
(with Eo fixed to 670 MeV). In order to test for curva- 
ture in the 7 -ray spectrum of PKS 2123—463 we used as 
an alternative spectral m odel the log-parabo l a (LP), dN/dE 
oc (E /Eq)~ qi ~ /31 ° s ^ e ^ Eo ' ) (lLandau et al.l[l986l : iMassaro et al.l 
l2004h . where the parameter a is the spectral slope at the en- 
ergy Eo and the parameter /3 measures the curvature around 
the peak. We fixed the reference energy Eo to 300 MeV. The 
fit with a log-parabola results in a TS = 554, with spectral 
parameters a = 2.35 ± 0.09, and /3 = 0.06 ± 0.04. Applying a 
likelihood ratio test to check the PL model (null hypothesis) 
against the LP model (alternative hypothesis), the power- 
law spectral model is favoured, indicating that no significant 
curvature was observed in the average 7 -ray spectrum. 


4 F. D’Ammando, A. Rau, P. Schady, et al. 


Figure [l] shows the 7-ray light curve of the first 41 
months of Fermi observations of PKS 2123—463 built us- 
ing 2- month time bins, with the exception of the final (1- 
month) data point. For each time bin the photon index was 
frozen to the value resulting from the likelihood analysis 
over the entire period. If TS < 10 2 -ct upper limits were 
calculated instead. The systematic uncertainty in the flux is 
energy dependent: it amounts to 10% at 100 MeV, decreas- 
ing to 5% at 560 MeV, a nd increasing to 10% above 10 GeV 
llAckermann et ahll2012ll . 

As shown in Fig. [TjPKS 2123—463 was in a low 7-ray 
state (0.1-100 GeV flux < 5xl0 -8 ph cm -2 s -1 ) for the 
first 2.5 years. In particular the source was not detected 
with 2-month time bin during the period 2009 August- 
2010 May. A first increase of the 7-ray flux was observed 
in 2011 February-March, and subsequently a flaring activ- 
ity in 2011 December. We show a light curve focused on 
the period of high activity (2011 December 4-2012 Jan- 
uary 4; MJD 55899-55930) with 2-day or 4-day time bins 
(Fig. [2j) . The peak of the emission was observed between 
December 13 15:43 UT and December 14 15:43 UT, with 
a flux of (128 ± 23) xl0~ 8 photons cm -2 s^ 1 in the 0.1- 
100 GeV energy range, a factor of ~ 25 higher with respect 
to the average 7-ray flux observed during 2008-2011. The 
corresponding observed isotropic 7-ray luminosity peak in 
the 0.1-100 GeV energy range is 8.9±0.8 x 10 48 erg s _1 
(assuming z = 1.46±0.05, see Sect. [SJ), comparable to the 
values r eached by the most powerful FSRQs in flari ng activ- 
ity fe.g. lOrienti et ahll2012l : lAckermann et aI1l2010l b Leav- 
ing the photon index free to vary during the period of high 
activity and Eq fixed to 670 MeV the fit results in a pho- 
ton index Y = 2.26 ± 0.06, showing a moderate harder- 
when- bri ghte r behaviour already observed in other FSRQs 
dAbdo et al.ll2010bj ). Replacing in the same period the PL 
with a LP, fixing the reference energy Eo to 300 MeV, we 
obtain spectral parameters a = 2.09 ± 0.13, and /3 = 0.10 ± 
0.04. We used a likelihood ratio test to check the PL model 
(null hypothesis) against the LP model (alterna tive hypoth- 
esis). These values may be compared, following iNolan et al.l 
d2012h . by evaluating the curvature Test Statistic "I'S „ 
= 2(log(like)i,p-log(like)pi J )=3 (~1.7-cr); thus no significant 
curvature was observed in the 7-ray spectrum also during 
the high activity period. 

Finally, using the gtsrcprob tool, which calculates the 
probabilities an event belongs to each of the sources in the 
model, we estimated that the highest-energy photon emitted 
by PKS 2123—463 was observed with a probability of 87% at 
distance of 0.35° from the source with an energy of 92 GeV. 
For a single high-energy photon (^ 1 GeV), sky location can 
only be determined within ~ 0.8° (95% confidence radius). 


3 SWIFT DATA: ANALYSIS AND RESULTS 

The Swift satellite dGehrels et al.|[2004T ) performed two ob- 
servations of PKS 2123—463 in 2011 December triggered 
by the 7-ray flare detected by Fermi- LAT. As a compari- 
son we also analysed two additional observations performed 
during 2010-2011. The observations were performed with 
all three on -b oard i nstruments: the X-ray Telescope (XRT; 
iBurrows et al.ll2005l. 0.2-10.0 key), th e UltraViolet Optical 
Telescope (UVOT; Roming et al.ll2005l . 170-600 nm) and the 


Burst Alert Telescope fBAT; iBarthelmv et alJ[2005l . 15-150 
keV). 

The hard X-ray flux of this source is below the sen- 
sitivity of the BAT instrument for the short exposure 
time of these observations. Moreover, the source was not 
present in the Swift BAT 58-month hard X-ray catalog 


catalog qCusumano et al 


Baumgartner et al.1 2010) and the 54- month Palermo BAT 


201011 . 


3.1 Swift /XRT 

The XRT data were processed with standard procedures 
(xrtpipeline vO.12.6), filtering,and screening criteria by 
using the HEAsof t package (v6.11]j. The data were collected 
in photon counting mode in all observations, and only XRT 
event grades 0-12 were selected. The source count rate was 
low (< 0.5 cnt s _1 ); thus pile-up correction was not required. 
Source events were extracted from a circular region with a 
radius of 20 pixels (1 pixel ~ 2 , . , 36), while background events 
were extracted from a circular region with a radius of 50 pix- 
els away from the source region. Ancillary response files were 
generated with xrtmkarf, and account for different extrac- 
tion regions, vignetting and point spread function (PSF) cor- 
rections. We used the spectral redistribution matrices v013 
in the Calibration database maintained by HEASARC. 

Considering the low number of photons collected (< 
200 counts) the spectra were rebinned with a minimu m 
of 1 count per bin and the Cash statistic dCash I Il979l l 
was applied. We fitted the spectrum with an absorbed 
power-law using th e photoelectric absorption model tbabs 
ll Wilms et al JlOOOI 'I. with a neutral h ydrogen column fixe d 
to its Galactic value (2.34 xlO 20 cm~ 2 ; lKalberla et al.ll2005l b 
The fit results are reported in Table |T| and Fig. [4] 

A marginal increase (1.5 cr) of the X-ray flux with re- 
spect to the previous Swift/XRT observations was observed 
on 2011 December 15, soon after the 7-ray flare. The source 
remained at a similar flux level in X-rays on 2011 December 
19. 


3.2 Swift/UVOT 

UVOT UV/optical imaging was obtained during all four 
Swift observations of PKS 2123—463. The photometry was 
carried out on pipeline processed sky images downloaded 
from the S wift data centefl following the standard UVOT 
procedure (IPoole et al.ll2008l l. Source photometric measure- 
ments were extracted from the UVOT imaging data using 
the tool UVOTMAGHIST (vl.l) with a circular source extrac- 
tion region that ranged from S'/S-S^radius to maximize the 
signal-to-noise. In order to remain compatible with the effec- 
tive area cali b rations, whic h are based on 5”aperture pho- 
tometry (IPoole et al.ll2008lh an aperture correction was ap- 
plied where necessary. This correction was at maximum 5— 
6 % of the flux, depending on the filter. The UVOT pho- 
tometry is presented in Table 0 and Fig. []] Contempora- 
neous UVOT observations in 2011 December 15 found PKS 
2123—463 about 0.5 mag brighter in v-, u-, and wl-band, 0.7 


3 http: / /heasarc. nasa.gov/lheasoft/ 

4 http: //www. swift. ac.uk/swift_portal 
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Table 1. Log and fitting results of Swift/yL RT observations of PKS 2123—463 using a power-law model with Ah fixed to Galactic 
absorption. 


MJD 

UT Date 

Exposure Time 
[sec] 

r 

Flux (0.3-10 keV)“ 
[xlO -12 erg cm -2 s -1 ] 

55467 

2010-09-28 

2712 

1.47±0.30 

1.31T0.46 

55745 

2011-07-03 

2912 

1.61T0.28 

1.20T0.32 

55910 

2011-12-15 

4947 

1.68±0.18 

1.90T0.34 

55914 

2011-12-19 

5521 

1.26±0.18 

1.83T0.27 


Observed flux. 


Table 2. UVOT Photometry of PKS 2123-463. 


UT Date 



AB Magnitude 0 




uvw2 

uvm2 

uvwl 

u 

b 

V 

2010-09-28 18:29-23:30 

19.78 ±0.09 

19.34 ±0.10 

19.17 ±0.10 

18.66 ± 0.09 

18.61 ± 0.14 

18.59 ± 0.20 

2011-07-03 17:15-23:29 

19.83 ±0.20 

19.28 ±0.14 

18.75 ± 0.11 

18.47 ±0.12 

18.81 ± 0.25 

18.21 ± 0.26 

2011-12-15 20:45-01:25 

18.97 ±0.06 

18.58 ±0.06 

18.24 ±0.05 

17.95 ± 0.06 

- 

17.64 ±0.06 

2011-12-19 04:20-04:39 

19.17 ±0.12 

18.80 ±0.12 

18.51 ± 0.09 

18.04 ±0.07 

17.77 ±0.07 

17.82 ± 0.15 


Corrected for Galactic foreground reddening of Eb— v = 0.03 mag llSchlegel et al.ll) 99<*j ). 


mag in m2-band, and about 0.8 mag in w2-band compared 
to the UVOT observation performed on 2011 July 3. 

4 GROND DATA 

On 2011 December 18, at 01:21 UT, PKS 2123—463 was ob- 
served with th e Gamma-r ay Optical/Near-Infrared Detector 
f GROND: iGreiner et alj|2008h mounted on the MPG/ESO 
2.2-m telescope at La Si lla, Chile. Preliminary results have 
been reported in iRau et al.l (1201 ll J . GROND is a 7-channel 
imager that observes in four optical and three near-IR 
channels simultaneously. The data were reduced and anal- 
ysed wit h the stand ard tools and methods described in 
[Kriihler et al.l d2008T ). Calibration was performed against 
an SDSS stan dard star fie ld ( c/'r'i' z') and against selected 
2MASS stars (ISkrutskie et al.ll2006l ) ( JHKs ). This resulted 
in l-(j accuracies of 0.05 mag ( g'r'i'z ' , J, H), and 0.07 mag 
(A's). All magnitudes are correcte d for Galactic foregr ound 
extinction of Eb-v = 0.030 mag (ISchlegel et al.lll998l ) and 
are summarised in Table [3] 


5 RADIO DATA 

PKS 2123—463 is monitored by the TANAMI program 
(Tracking Active gal actic Nuclei with Austral Milliarcsec- 
ond Interferometry; lOiha et~ah IM) at a number of ra- 
dio frequencies and resolutions with the Australia Telescope 
Compact Array (ATCA) and the Ceduna facilities. ATCA 
is observing this source every few weeks with “snapshot” 
observations at frequencies from 5.5 GHz through 40 GHz 
where each frequency is the centre of a 2 GHz wide band 
and the fluxes are calibrated against the A TCA primary flux 
calibrator PKS 1934—638 (iStevens et al.if2012l ). Each flux 
density has a 1 -a uncertainty of 0.01 Jy. The Ceduna radio 
telescope in South Australia is monitoring PKS 2123—463 at 
6.7 GHz. Each flux den sity has a 1-a uncertainty of 0.3 Jy 
llMcCulloch et af1l2005l h No sign of increased activity of the 


Table 4. ATCA (5.5, 9, 17, 19, 38, and 40 GHz) and Ceduna (6.7 
GHz) observations of PKS 2123—463 in 2011. 


UT Date 

Frequency (GHz) 

Flux (Jy) 

2011-05-17 

5.5 

0.86 


9.0 

0.77 


17.0 

0.67 


19.0 

0.64 


38.0 

0.45 


40.0 

0.44 

2011-10-14 

5.5 

0.86 


9.0 

0.76 


17.0 

0.60 


19.0 

0.56 


38.0 

0.42 


40.0 

0.41 

2011-11-08 

5.5 

0.85 


9.0 

0.71 


17.0 

0.57 


19.0 

0.53 


38.0 

0.43 


40.0 

0.43 

2011-12-19 

6.7 

0.81 

2012-02-18 

6.7 

0.78 


flux density was detected between 40 and 5.5 GHz before 
the 7-ray flare (see Fig. 0 . 

In addition observations of PKS 2123—463 were made 
after the 7-ray outburst using the KAT-7 array (the proto- 
type for MeerKAT) in the Karoo. Since the array is still un- 
dergoing commissioning tests some uncertainties remained 
about the absolute calibration scale. To minimize these un- 
certainties the observations were taken at similar times of 
days (hence Local Sidereal Time range) in late 2011 Decem- 
ber to early 2012 January, with 4-5 hr durations. The very 
first observation (December 21) failed because the online 
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Table 3. GROND Photometry of PKS 2123—463. 


UT Date 

9' 

AB Magnitude 0 

r' i' z' 

J 

H 

K s 

2011-12-18 01:21 

17.78 ±0.05 

17.52 ±0.05 17.43 ±0.05 17.24 ± 0.05 

16.75 ± 0.06 

16.53 ± 0.06 

16.24 ±0.09 


Corrected for Galactic foreground reddening of Eb— v = 0.03 mag iSchlegel et al.lfl998h . 


fringe stopping was not working properly, but as the array 
is very small (longest baseline 200-m) and the frequency low 
enough (1.822 GHz) the delay tracking corrections could be 
done offline with only marginal loss in signal to noise ratio. 
Subsequent observations were done in this mode. There were 
also different problems with various receivers over this time 
so «t>coverage and maps were not directly comparable. The 
bandpass and gain solutions (on PKS 1934—638) were good 
enough and checks against a nearby phase reference source 
(PKS 2134—470) gave sufficient confidence in the visibilities 
and phases to make simple maps of the inner 0.5° of the 
primary beam and verify that no sources stronger than 100 
mjy were present. 

Assuming a flux density for PKS 1934—638 of 13.6 Jy 
at 1.822 GHz (as interpolated from ATCA models) the mea- 
sured flux density of PKS 2123—463 is 1.13 Jy on 2011 De- 
cember 24, 26, 29, and 2012 January 7 and 13. The uncer- 
tainties on the absolute flux density scale is 5%. As a com- 
parison we extrapolate the flux density at 1.82 GHz using 
the values at 2.7, 0.8, and 0.4 GHz from the Parkes Cata- 
logue (1990) and the Molonglo Sky Survey (2008). Assuming 
the spectral index 0.4 derived between 2.7 and 0.4 GHz, we 
obtain a flux density at 1.82 GHz of 1.0 Jy, indicating that 
no significant variation was observed during the 7-ray flaring 
activity at this frequency. 

6 NEW PHOTOMETRIC REDSHIFT 

ESTIMATION 

The redshift of PKS 2123—463 has not been established con- 
vincingly yet. Following the method described in iRau et akl 
(l2012ll and applied to a sample of 103 Fermi - LAT blazars, 
we combined the GROND photometry with the UVOT pho- 
tometry from 2011 December 19 in order to construct a 13- 
band spectral energy distribution (SED). In order to ac- 
count for source-intrinsic variability between the GROND 
and UVOT observing epochs, the GROND photometry was 
corrected by —0.05 mag in all bands fsee IRau et alj I2OI2L 
for a detailed description of the method). The resulting 
SED, corrected for the Galactic foreground reddening of 
Eb-v = 0.03 mag l|Schlegel et al.lll998 ). has been fit with 
a set of power-law models as well as with hybrid tem- 
plates built from normal galaxies and AGN (ISalvato et al.l 
2009|j201l|) using the Le PHARE code (lArnouts et al.ll 19991 : 
Ilbert et al l 12006! ) . The best-fitting redshifts for both tem- 
plate libraries are in good agreement (see Fig. [3]). For the 
power-law model we obtain a 99% confidence redshift of 
a = 1.37 !q' 2 i (x 2 — 17, Pz = 94.6 0) and for the hybrid 

5 The best fit templates and redshifts were selected on the ba- 
sis of a simple y 2 method. P z is the integral of the probability 
distribution function j f(z)d,z at z p hot ± 0.1(1 + z p hot)i which 



Wavelength (pm) 

Figure 3. UV-near-IR SED composed of Swift / UVOT observa- 
tions obtained on 2011 December 19 and GROND data taken on 
2011 December 18. The simple power-law model (solid line) and 
the AGN /galaxy hybrid template (dashed line) suggest a photo- 
metric redshift of z ^ 1.45 (see text for details). 

models we find z = 1.46 ± 0.05 (y 2 = 23, P z = 99.2). This 
is in rough agree ment with the init i al spe ctroscopic redshift 
of z = 1.67 from lSavage fe Wright! (ll98llJ . 


7 SED MODELING AND CONCLUSIONS 

Beyond the excellent spatial association obtained with the 
41-month Eermi-LAT data set, the most secure and dis- 
tinctive signature for firm identification of the 7-ray source 
detected by Eerrm-LAT with the blazar PKS 2123—463 is 
the simultaneous increase observed in the 7-ray, X-ray, and 
optical-UV bands (see Fig. [4] Sections [2] and [3]). 

In order to investigate the physical properties of the 
source we have built a simultaneous SED of the flaring state 
of PKS 2123—463. The Eermi-LAT spectrum was built with 
data from observations centred on 2011 December 10 to 19 
(MJD 55905-55916). In addition we included in the SED 
the GROND and Swift (UVOT and XRT) data collected 
on 2011 December 18 and December 19, respectively. Here, 
the flux suppression, in particular in the UV bands, was 
corrected assuming a power-law spectral shape consistent 
with the optical-near-IR measurements (F\ oc A 0,8 ). The 
data from ATCA and KAT-7 collected on December 19 and 
24, respectively, provided information about the radio part 
of the spectrum. The flare centred on ~ MJD 55909 had 
a variability timescale of ~ 2 days, which constrains the 
size of the emitting region during the flare to Rj, <2.1 x 
10 16 (<7d/20) cm. 


describes the probability that the redshift of a source is within 
0.1(1 + z ) of the best fit value. 
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Figure 4. Light curves of PKS 2123—463 collected in 7 rays by 
Fermi - LAT (0.1-100 GeV flux in units of 10 -8 ph cm -2 s _1 ), in 
X-rays by Swift /X RT (flux in units of 10 12 erg cm -2 s _1 ), in UV 
(filters wl: empty squares, m2: empty triangles, w2: filled circles; 
flux in units of 10 12 erg cm -2 s — 1 ), and in optical (filters u: 
empty squares, b: empty triangles, v: filled circles; flux in units of 
10 13 erg cm -2 s -1 ) by Swift/ UVOT, and radio by ATCA (empty 
squares: 40 GHz, empty triangles: 17 GHz, filled circle: 5.5 GHz) 
between 2010 June and 2011 December. 


torus is the likely seed photon source (e.g., Marscher et alj 
l2010h . From the disc luminosity obtained by the fit of the 
blue bump, Ldisc = 1-8 x 10 46 erg s _1 , we estimate an asso- 
ciated BLR radius Rblr = 4.4 x 10 17 cm, based on the re- 
lation between disc luminosity and Rblr determined from 
reverberation mappin g campaigns (e.g.. iKaspi et all 120051 : 
iGhisellini et al.ll2008lh To minimize the scattered BLR con- 
tribution, we placed the emitting region at r > Rblr ■ Here 
the primary seed photon source is the dust torus, which 
was simulated as a one-dimensional ring with radius Rdust 
aligned orthogonal to the jet, emitting as a blackbody with 
temperature Tdust and luminosity Ldust- 

The model fit to the broadband SED can be seen in 
Fig.[5]and the parameters can be found i n Ta ble[5] A descrip- 
tion of the parameters can be found in iDernrer et ali (120091) . 
The dust parameters were chosen so that R dust is roughly 
consi stent with the sublimation radius (iNenkova et ali 
2008). The EC-BLR component was calculated assuming the 
seed photons are from Ha and have a luminosity of 0.1 Ldisc- 
The electron distribution used, a broken power-law with 
index pi = 2.0 below the break at j' brk and p 2 = 3.8 for 
Ibrk < : is consistent with particles injected with index 

2.8, and emission taking p lace in the fast cooling regime (e.g., 
iBottcher fe Dermeifcool h That is, particles are injected be- 
tween and 7 ' max , with a cooling electron Lorentz factor 

3 m e c 2 

Icool = - 7 

4 CdTiescUtot 


A “blue bump” accretion disc component is clearly vis- 
ible in the optical/UV data. We modeled these data with 
a combination of a nonthermal synchr otron compone nt and 
a Sh akura-Sunyaev disc component (IShakura and Sunvaev I 
1 1973T) . A fit to the disc component allows us to get a rough 
estimate of the black hole mass of Mbh ~ 2x 10 9 Mq. This is 
the first black hole mass estimate for this source. It is con- 
sistent with black hole estima tes for o ther liig h-2: F SRQs, 
obtained in a similar way by IGhisellini et al.l (|201l|). Our 
mass estimate for PKS 2123—463 follows Ghisellini et al.l 
(120111) and assumes that the innermost stable circular or- 
bit (ISCO) is Rdisc,i sco = 6 R g , as one would expect for a 
Schwarzschild (i.e. , nonrotating) black hole. If the jet is pro - 
duced by the Blandf ord-Znaje k (IBlandford fc Znaiekiri977l l 
or Blandford-Payne (IBlandford fe Pavnel 1982h mechanisms 
this requires a nonzero black hole spin, and one expects 
Rdiscjsco t° be smaller or larger, depending on w heth er th e 
spin is retrograde or prograde (e.g., iGarofalo et al.l |2010 i) . 
Due to the uncertainty in spin, this mass estimate can be 
considered to have considerable uncertainty. 

We modelled the portion of the SED from X-rays to 7 
rays assuming emission from a relativistic jet with mecha- 
nisms of synchrotron self-Compton (SSC), and Compton- 
scattering of a dust torus external to the jet (EC-dust). 
The d escri ption of the m o del c an be found in iFinke et al.l 
ll2008f) and iDermer et al.1 (120091 4 . The synchrotron compo- 
nent considered is self-absorbed below ~ 10 11 Hz. Correla- 
tions of 7-ray and optical flares with radio light curves and 
rotations of optical polarisation angles in low-synchrotron- 
peaked blazars seem to indicate the 7-ray/optical emitting 
region is outside the broad line region (BLR), where the dust 


where utot is the total energy density in the blob frame, 
which in the case of our model fit is dominated by the ex- 
ternal energy density. I11 this case 7 coo i is associated with 
7 min, since it is in the fast cooling regime. Also note that in 
this model fit the magnetic held and electrons are nearly in 
equipartition. Jet powers were calculated assuming a two- 
sided jet. 

The Compton dominance for PKS 2123—463, i.e., 
the ratio of the peak luminosities of the Compton and 
synchrotron components, is « 50, which is a ra ther 

stan dard value for powerful blazars (e.g., IGhisellini et ali 
[ml]). A large disc luminosity was estimated from the 
UVOT data, as expected for a powerful FSRQ, with a 
Ldisc /IsE dd =0.2 in a greement w ith the blazar divide pro- 
posed by IGhisellini et ahl (j2009l f as a result of the changing 
of the accretion mode. 

Variability is common in 7-ray blazars and provides 
a powerful tool to associate them definitively with objects 
known at other wavelengths and to study the emission mech- 
anisms at work. The combination of deep and fairly uniform 
exposure over ~3 hr, very good angular resolution, and sta- 
ble response of the Fermi - LAT is producing the most sen- 
sitive, best-resolved survey of the 7-ray sky. On the other 
hand, those cases where there is a decrease in the activity, 
and thus of the significance of detection, can lead to a more 
complex identification process for a 7-ray source. When com- 
bined with simultaneous ground- and space-based multifre- 
quency observations, the Fermi-LAT achieves its full capa- 
bility for the identification of the 7-ray sources with coun- 
terparts at lower energies and the knowledge of their emis- 
sion processes, as reported here for the high-2, Compton- 
dominated FSRQ PKS 2123—463. 
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Table 5. Model parameters for the SED shown in Fig. [5] 


Redshift 

Bulk Lorentz Factor 
Doppler Factor 
Magnetic Field 
Variability Timescale 
Comoving radius of blob 
Jet Height 

Low-Energy Electron Spectral Index 

High-Energy Electron Spectral Index 

Minimum Electron Lorentz Factor 

Break Electron Lorentz Factor 

Maximum Electron Lorentz Factor 

Disc Luminosity 

Black Hole Mass 

Accretion efficiency 

Gravitational Radius 

Inner Disc Radius 

Outer Disc Radius 

Dust Torus luminosity 

Dust Torus temperature 

Dust Torus radius 

Jet Power in Magnetic Field 

Jet Power in Electrons 


z 

1.46 

r 

20 

<^D 

20 

B 

0.8 G 

ty 

1.7xl0 5 s 

K 

4. lx 10 16 cm 

r 

l.OxlO 18 cm 

Pi 

2.0 

P2 

3.8 

y ■ 

'ruin 

4.0 

^ brk 

1.0x10 s 

"7 max 

1.0x10 s 

L disc 

1.8xl0 46 erg s - 

M bh 

2 x 10 9 M 0 

V 

1/12 

R 9 

1.76X10 14 cm 

R disc.ISCO 

6 Rg 

R disc-max 

Oi 

o 

Ldust 

1.5xl0 45 erg s 

Tdust 

1.7xl0 3 K 

Rdust 

3.2xl0 18 cm 

Pj,B 

3.3X10 45 erg s“ 

p. 

1 j,par 

1.5xl0 45 erg s _ 



Figure 5. Spectral energy distribution data (circles and squares) 
and model fit (solid curve) of PKS 2123—463 with the model 
components shown as dashed curves. The data points were col- 
lected by GROND (2011 December 18), Swift (UVOT and XRT; 
2011 December 19), and Fermi-LAT (2011 December 10-19), to- 
gether with radio data from ATCA (2011 December 19) and KAT- 
7 (2011 December 24). 
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